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Abstract 

o ■ 

We investigate the off-forward matrix element of the light cone vector oper- 
ator for a dressed quark state in light-front Hamiltonian perturbation theory. 
We obtain the corresponding splitting functions in a straightforward way. We 
i—j | show that the end point singularity is canceled by the contribution from the 

qq ' normalization of state. Considering mixing with the gluon operator, we verify 

the helicity sum rule in perturbation theory. We show that the quark mass 
. effects are suppressed in the plus component of the matrix element but in the 

transverse component, they are not suppressed. We emphasize that this is a 
| particularity of the off-forward matrix element and is absent in the forward 

7— I 1 case. 

: 

o 
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^ ' 1. Introduction 

a: 

The hadronic matrix elements of quark and gluon operators appear in the description 
• of all scattering processes. They are in general of two types: in the inclusive processes, one 
>- ! encounters diagonal matrix elements of bilocal operators. These matrix elements are related 
to parton distributions. On the other hand, in the elastic exclusive processes, one encounters 
form factors, which are off-diagonal matrix elements of local operators. The generalized 
parton distributions interpolate between these two types of matrix elements [1]]. These are 
off-diagonal matrix elements of light-front bilocal operators. They play an important role 
in the deeply virtual Compton scattering amplitude || and electroproduction of mesons 
(for reviews of hard exclusive reactions, see ||). The off-forward matrix elements are 
the generalizations of the above two types of matrix elements; parton distributions are the 
forward limits of generalized parton distributions (GPD) and form factors are moments of 
them. 

Recently, the generalized parton distributions have been investigated in the light-front 
formalism by several authors and an overlap representation for the plus component in terms 
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of light-front wave functions has been given |7]||. Also GPD's have been constructed us- 
ing light cone model wave functions |J. The perpendicular and the minus components are 
somewhat more complicated because the operators in these cases involve interactions. They 
are the higher twist components. In this work, we calculate the the off-forward matrix ele- 
ments of the plus component and the mass dependent helicity-flip part of the perpendicular 
component of the bilocal vector operator for a dressed quark target in light-front Hamilto- 
nian QCD. Recently, we have made considerable progress in understanding polarized and 
unpolarized deep inelastic scattering (DIS) structure functions in this approach and we 
have shown that it gives an intuitive picture of DIS. It is suitable to calculate the forward 
matrix elements of transverse and minus components of the bilocal current operator. Also, 
the presence of quark mass does not cause any problem. Interference effects are straight- 
forward to handle. The splitting functions are obtained easily and they agree with the well 
known expressions jlOfl . It is possible to derive new sum rules which connect DIS structure 
functions to light-front QCD Poincare generators. In this work, we extend our approach to 
off-forward matrix elements. 

2. Plus Component 

We work in the so called symmetric frame |7y§[]. The momentum of the initial state is 
P^ L and that of the final state is P' M . The average momentum between initial and final state 
is then P» = ^±^. 

The momentum transfer is given by = P' M — P^, P' ± = —P± = 4r, skewedness 

£ = — ^qr. Without any loss of generality, we take £ > 0. We also get A~ = ^r. 
The off-forward matrix element is given by, 




(a) (b) 

Fig. 1: Light-front time ordered diagrams considered in the calculations of Eqs. ([!]) and 
(E9h; only one time ordering is important in the kinematical region considered. 
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The + component of the matrix element is parametrized in terms of the off-forward dis- 
tributions H(x,£,t) and E(x,£,t). However, the matrix element can be expressed directly 
in terms of overlaps of light-front wave functions. We calculate the plus component of this 
matrix element taking the target to be a dressed quark to order g 2 (see Fig. 1 (a)). We work 
in the light-front gauge, A + = 0, where the path-ordered exponential between the fermion 
fields in the bilocal operator is unity. For simplicity we suppress the flavor indices. 
The Fock space expansion of the operator is given by, 

. „ r dk + d 2 k x f dk' + d 2 k' L 

O = 4 > / t: / rr 

* J 2(2ir) 3 V¥ J 2(2nyVk 

[5{2xP + - k' + - k + p(k, s)b{k', s) 

+ 5{2xP + + k' + + k + )d{k, -s)d\k', -s) 
+ 5(2xP + + k + - k' + )d(k, -s)b(k', s) 

+ 5{2xP + + k' + - k + )tf(k, s)d\k', -s)] . (2) 

We have, k + > 0,k' + > 0, k + — k' + = p + —p' + = 2£ > p + . In the kinematical region, £ < x < 1, 
only the first term in Eq. (0) contributes . We restrict ourselves to this kinematical region. 

We take the state | P, a) to be a dressed quark consisting of bare states of a quark and 
a quark plus a gluon: 

| P,a)=cf h tf(P,a) | 0) 



G\ ,A 



+ E 7 ' - 2 \:\ J2(2^P+5\P -h-k 
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p{2Tlfk+ J V /2(27T)3fc 2 4 



<j> 2 (P,cr | k 1 ,a 1 ]k 2 ,X2)b t (k 1 ,a 1 )a^(k2,X 2 ) \ 0). (3) 

Here a' and are bare gluon and quark creation operators respectively and 4>i an d 4>2 are 
the multiparton wave functions. They are the probability amplitudes to find one bare quark 
and one quark plus gluon inside the dressed quark state respectively. Up to one loop, if one 
considers all kinematical regions, there will be non-vanishing contributions from the overlap 
of 3-particle and one particle sectors of the state, this situation is similar to QED |7|. In the 
kinematical region we are considering, such kind of overlaps are absent and it is sufficient to 
consider dressing only by a single gluon. The state is normalized to one. 0i actually gives 
the normalization constant of the state |]10[: 



1 2 



>1 



a s „ f 1 ^ 1 + x 2 Q 2 



1 - -±C f / dx— log^j, (4) 

2f ie 1 - X U 1 



fi 

within order a s . Here e is a small cutoff on x. 
The matrix element becomes, 

F+ = A\fP + Vp' + [<t)l{P')4>i{P)5{2xP + - 2P + ) 

+ J2 J dptd 2 pi ( f ) ;(P , 1 A+p 1 ,P-p 1 ) ( j )2 (P,p 1 ,P-p 1 )5(2xP + -A + -2pi)] (5) 

X] denotes summation over helicities of the quark and gluon. The first term comes from one 
particle sector and it contributes only when x = 1. The first term receives contribution upto 
order a, from the normalization condition of the state. 
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The other nontrivial contribution comes from the two particle sector given by the second 
term. 

We introduce Jacobi momenta Xj,^ -1 such that J2i x i — 1 an d J2iQi' L — 0. Also, we 
introduce boost invariant wave functions, 



fa = 0!, M*i,<lt) = VP (/>(kt,ki x ). (6) 

The contribution from the two-particle sector then becomes, 

F+ = 2^/dV^(f^|,g ± + ^|A^ 2 (f±|,^). (7) 

The two particle wave function depends on the helicities of the quark and gluon and is given 
in terms of tpx as, 

in x(l — x) ^ n 1 q x r q 1 - a 1 - ■ q 1 - , 



'1 — X 



ima~ 



x 



xA*i>i- (8) 



Here a 1 = a 2 and cr 2 = —a 1 . Using Eq. (|7|) and Eq. (|8[), we see that the mass terms give 
suppressed contributions. Since mass terms in the vertex cause helicity flip, helicity flip 
parts of the matrix element are suppressed. We calculate the helicity non-flip part. 
Using Eq. (|j) we get, from Eq. (^), 

W?' C ' (1-x) ( 9W + ^>A^ (1 - 2 « (9) 

where C/ = ^jv^ for SU(N). The g 1 - integral is nontrivial and it is divergent for large q^. 
Integration over the polar angle gives, 

rf 9 7 i i a = 2tt / gdg 



(g ± ) 2 (g ± + aA ± ) 2 7 M ^ q 2 + a 2 A 2) 2 _ 4a 2 A 2 ?2 

A r (<? 2 + a 2 A 2 ) 



7T 



A* 9 



( ? 2 + fl 2 A 2)2 _ 4a 2 A 2 g 2 



- 1 



= h + h, (10) 

where a = rfEpy- A * s the large transverse momentum cutoff and \l is the factorization scale 
separating hard and soft dynamics ||10|| . The divergence structure of the above integral can 
be seen by expanding the denominator of the integrand in the limit of small A =| A^ |. We 
find that 1\ is logarithmically divergent in the limit A — > oo but Ii has no divergent part. 
The full divergent part is given by, 

P+ 9 9 2 r v^?(l-2g 2 + x 2 ) Q 2 
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where we have cut the transverse momentum integral at some scale Q 2 . The splitting 
function can be easily extracted from the above expression : 



(\ + - 2f 2 ) 

^■fl = c ' ( i-s)(i-V) - (12) 

This agrees with || (with | replaced by £). The above expression contains end point 
singularity at x = 1. However, this is canceled by the contribution from the normalization 
of the state to the single particle matrix element, so that the final result becomes [], 

F+ = - *> + T, C '^% S{1 - ~ x) + (13) 

This result shows the importance of the normalization contribution to the single particle 
matrix element. In other words, it includes contributions from virtual gluon emission. In 
this approach, one uses probability amplitudes rather than probability densities in Altarelli- 
Parisi method and the effects due to both real and virtual gluon emissions are taken into 
account to the same order in a s without any difficulty In order to get the full scale evolution, 
one has to consider all the kinematical sectors which is beyond the scope of the present work. 
Also, one can see that the final result has no singularity at x = £ which is as expected [[?]. 

Next, we calculate the helicity flip part of the matrix element. The helicity flip contri- 
butions come from the mass term in the expression of the two-particle wave function. The 
form of the wave function shows that this contribution is suppressed. 

Next, we parametrize the off-forward matrix element in terms of the generalized quark 
distributions, 



1 



% 



F^ = W U y (P') H q (x,U)r + E q (x,Z,t)—a» a A a U X (P) + (14) 



2M 



where the ellipses indicate higher twist terms. U\(P) is the quark spinor in our case. Using 
the explicit form of the light front spinors, we get for the plus component, 
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F + S x , x = 2y/l-eH q (x,£,t) - ^L= E q (x,C,t). (15) 

Here we have calculated the helicity non-flip part. 
The helicity flip part becomes, 

-A 1 + iA 2 

F+5 X ,- X = ^—. Eq {x^t). (16) 

mVl — ? 

Since the helicity flip part of the matrix element is suppressed, we naturally find that E is 
suppressed in perturbation theory. Here we have taken both m and | A 1 - | to be small. In 
the limit of small £, we expand H q (x,£,t) ~ H q (x, 0, 0) + £, 2 H' q (x, 0, 0) in the kinematical 



Here (i^ x ) + is the usual (principal value) plus prescription. 
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region £ < x < 1. Using Eq. (|i"3|) and Eq. ([15]) and equating the coefficients of equal powers 
of £ on both sides, we get, 

rv <9 2 /S (1 4- f 2 ^ \ 

H q (x, 0, 0) = 5(1 - x) + ^C f log^ (±6(1 - x) + (17) 



H' g (x,0,0) = --^C f log^(l+x). (18) 



Eq. ( |iT| ) gives the forward limit of the generalized quark distribution H q and it can be easily 
identified with the unpolarized quark distribution for a dressed quark state in perturbation 
theory p[. 



Subsequently, we calculate the gluon distribution, 

= ~~ J ^^(P'X' | F+»{- Z -)F+{^) | PA). (19) 
The Fock space expansion of the relevant part of the operator is given by, 

O g = _| + (2(2 ^ )3)2 E / dk+d 2 ki J dktd 2 k^(k 1 ,X)a(k 2 ,X)5(2xP + - kf - k+) (20) 

We calculate the matrix element for a quark state dressed with a gluon (see Fig. 1(b)). The 
Fock space expansion of the state is given by Eq. @. 
The matrix element is given by, 

F 9=\T. J WK^riM^,* 1 - + T^A^V^e (21) 

Using the full form of the two particle wave function, we find that the helicity flip terms 
proportional to the quark mass give suppressed contribution and the helicity non-flip part 
is given by, 

:a Q 2 (1 + (1 - x) 2 - C 2 ) 



K = 2 V 1 - ^C f log ^- ^' ^ . (22) 



The splitting function can be easily extracted and agrees with |2| with | replaced by £. 

The gluon matrix element is parametrized in terms of the twist two distributions, H g 
and Eg, 

F^x = ^Uy(P')[H g (x,tt)r + E g (x,^t)^a^A a ]u x (P) + (23) 

The fact that the helicity flip part of the matrix element is suppressed means that E g is also 
suppressed. In the limit of small £, we again expand H g (x, £, t) ~ H g (x, 0, 0) + £ 2 H' g (x, 0, 0) 
in the kinematical region £ < x < 1 and we obtain, 

H g (x, 0, 0) = ^C f log% {1 + {1 _~ m , (24) 
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One can identify the above expression with the unpolarized gluon distribution for a dressed 
quark target JlO . 



H>(xAO) = ^C f log%^-jy. (25) 

Z7T fl z X 



We next verify the helicity sum rule in perturbation theory. It has been shown JTTJ 



that the light front helicity operator when expressed entirely in terms of the dynamical 
fields in the light front gauge, has the same form as in the free theory, provided we restrict 
ourselves to the topologically trivial sector (i. e. we take the dynamical fields to vanish 
at the boundary). This eliminates the residual gauge degrees of freedom and removes the 
surface terms. The helicity operator is given by, 

J 3 = J% + J% + J 3 + J 3 go , (26) 

where, J 3 j is the intrinsic quark helicity, J 3 Q is the orbital quark helicity, J gi is the intrinsic 

f'' 

go 



gluon helicity and J 3 is the orbital gluon helicity. The operators are given by, 



J 3 fo = J dx-d 2 x L ^ +] i{x l d 2 - x 2 d 1 )iP + , (27) 
Jfi = -J dx~d 2 x ± ij + ^ 3 ^ + , (28) 
J 9 3 o = I J dx~d 2 x ± [x 1 {d + A 1 d 2 A 2 + d + A 2 d 2 A 2 ) - x 2 (d + A 1 d 1 A 1 + d + A 2 d 1 A 2 )], (29) 

= 1 J dx~d 2 x ± (A 1 d + A 2 - A 2 d + A 1 ). (30) 

The color indices are implicit in the above expressions. We can check explicitly using the 
above expressions that the helicity sum rule for a dressed quark in perturbation theory is 
given entirely in terms of H(x, 0, 0), 

Jo 



a O 2 A 2 2 

dxxH q (x, 0, 0) = 1 - -±C f log^- = -{P, fl J% \ P, T) + ^{P, T| J% \ P T), (31) 

z7r ji z 3 I\ Jy 



J Q dxxH g (x, 0, 0) = ^LCflogZp- = -(P, T| J gi \ P, T) + -^(P, T| J 3 go \ P T), (32) 

where iV is a normalization constant. This gives, for a dressed quark state in perturbation 
theory, 

J o 1 dxx(H q (x,0,0) + H g (x,0,0)) = 1(P,T| 2 J 3 | P, t) = 1. (33) 



By comparing the Ihs and rhs of Eqs. ([H]) and (p2|) , one also verifies that in perturbation 



theory the total quark (gluon) momentum contribution is identical to the total quark (gluon) 
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helicity contribution. This result is expected from the fixed point solutions of the leading 
log evolution equations of J q and J g calculated in [12], which shows that the partition of the 



nucleon spin between quarks and gluons follows the partition of nucleon momentum. 
3. Helicity Flip Part of the Transverse Component 

Having studied the plus component of the matrix element of the bilocal operator, we 
now show the importance of quark mass in the matrix element of the transverse component 
of the bilocal vector operator. The matrix element of the transverse component is given by: 

i& = / ^ P+ ^(P'X' | <A(-y)7>(y) I FA). (34) 
The bilocal operator in this case can be written as, 

O = <K- y hV( y) = ^ +t (- ( y ) + r f (- y)<^ + ( y)- (35) 

The operator involves the constrained field and therefore it is the so-called bad 

component. ip~ can be eliminated using the constraint equation in light-front gauge, 

id+i/j- = [a 1 - ■ (id 1 - + gA 1 ) + 7°m]^+. (36) 

The terms linear in mass produce helicity-flip contributions in the matrix element. Quadratic 
mass terms generate helicity non-flip terms but it can be shown that they are suppressed. 
It is known that in the forward case, for non-zero P -1 , the above matrix element is related 



to the unpolarized DIS structure function F 2 , as calculated through the plus component ]T3" 
The contribution of the mass dependent helicity-flip part of the operator in this case cancels 
between the two terms of Eq. (^), as a result, quark mass effects are still suppressed here 
in the forward limit. 

We calculate the matrix element of the helicity flip part of the operator in the off-forward 
case. 

The operator is given by, 

O = O m + O fc x + O g , (37) 

where O m is the explicit mass dependent part of the operator, O k ± is the explicit k± depen- 
dent part and O g gives the interaction dependent part. The Fock space expansion of O m is 
given by, 

o m = -^(-yb^+(^) + ["-^(-^h^) 



2E 



dk+cPk 1 r dk'+d 2 k 



2U± 



2(2tt)VF J 2(2vr) 3 \/£ 



(im)5(2xP + -k + - k r+ p(k, s)b(k', s ') x y Xs/ (^ ~ ^:)- (38) 

Here \ is the two component spinor. This is the part of the operator that is relevant in 
the kinematical region £ < x < 1. We take the state to be a quark dressed with a gluon as 
before. The matrix element is given by, 
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2£ 1 

= ~{im)xU 2 X(Tl -j=== 



+ d 2 q A 



1 



5(1 - x)ip*ii>i 
1 — x 



.^ + r ^A> 2 (f±|,/) 



(39) 



Using the explicit form of the two particle wave function, and also using the normalization 
condition of the state, we write this as, 



6(1 



°> • Q 2 1 3 <■ , — \ - 2£ 2 



2tt~ j "~ ; V V2~ v ~ " y ' (1 -x) + (x 2 -£ 2 )- 
The other contribution to the helicity flip matrix element comes from, 



a, - £t ( _£l )(5 i + ia,d 2 )^(\) + [(^ - ^3d 2 )^ f (- y 



^ = e 

Here £ is the two component fermion field. The Fock space expansion is given by, 
dk+cPk 1 - r dk'+d 2 ^ 



(40) 



(41) 



-/+ 



2(27r) 3 v^ J 2(2nfVk 
S(2xP + - k' + - k + )tf(k, s)b(k', s) — - i<r 3 e ij ') + — + *a 3 e ij ) 



(42) 



The two particle contribution to the matrix element is of the form, 



1 



2x 



+ T^E/ <*V«(f^f .f 1 + ^A^Wf^Jx^V)*^. (43) 



We have taken | A -1 | to be small. 

The terms linear in mass in ^2 give the helicity flip contribution given by, 



F k^ = -5z(i™>)X> Xa> 



-rC f log 



Q 2 a s 2(l-x)(x 2 + e + 2x) 



p+^—'™- ~" J ^ 2 2<K (x 2 -e)(l~C 2 ) 

The interaction part of the operator is given by, 
Fa = 3t\- Z -7r)^r^ + ia,A 2 )i(^-) + g 



2 ' id ■ ' ' •' " v 2 
The Fock space expansion of the operator: 



-^(- T )(^-^( T ). 



(44) 



(45) 



Fi= E 



dk-^ d k^ f dk<2 d k^ f dk^ d k- 



Sl , S2 ,x J 2(2ir) 3 Vkl J 2(2tt) 



45(2xP + -kf-k+ + fc+) 



«0"3e 2 )x 



A' 2 



+ 45(2xP + - Jfe+ - jfc+ -k 



3 



*2 



+& t (/ci, si)6(fc 2 , ^y^s, A) 45(2rcP+ - kf - k+ - k 



^2 "I" 



*2 



3 



+ 4<5(2xP+ - fcf - fc 2 + + H f 1 ^, + 'TP* 32 } ) ■ (46) 



The terms containing d and d) will contribute in higher order. Out of the four terms, only 
the first and the fourth terms can contribute to the matrix element. These contributions 
are given in terms of the overlap of two-particle and one-particle wave functions. However, 
these contributions are zero since J2\xi( e \^- ~ i <J 3 e \)(.c rl ~ ' e \*))Xa' — 0- 

So we get, from Eq. (JHJ) and (f44|), the helicity flip part of the matrix element: 



F = - — {im)xl<7 X*> 



6{l-x) 



p+ k //vct ^ vr^p 

+ Y, c ^k 5{l - x) + 

- ¥^°V (i-e)^-e) J ■ (47) 

So we see that in contrast to the forward case, the effect of quark mass is not suppressed 
in the matrix element of the transverse component of the bilocal current. It is clear that 
such a helicity flip contribution is zero in the forward limit because the above expression is 
proportional to £. So, this effect is a particularity of the off-forward matrix element only. 
An interesting study using this approach will be to investigate the Wandzura-Wilczek 



relation |14[ for the off-forward matrix elements of the transverse vector and axial vector 
operators in perturbation theory, for dressed quark states. We plan to undertake such studies 
in the near future. 

To summarize, in this work, we have calculated the off-forward matrix element of the 
light-cone bilocal vector operator for a dressed quark state in perturbation theory. We have 
restricted ourselves in the kinematical region £ < x < 1. The contribution from the over- 
lap of three and one particle wave functions is absent in this case. We have obtained the 
corresponding splitting functions directly. The end point singularity is canceled by the con- 
tribution from the normalization condition of the state. We have shown that the generalized 
parton distributions E q and E g are suppressed in perturbation theory. Furthermore, we have 
verified the helicity sum sule in perturbation theory for a dressed quark state. The terms 
linear in quark mass cause helicity flip. However, such terms are suppressed in the matrix 
element of the plus component. We have calculated the helicity flip part of the matrix el- 
ement of the transverse component of the same operator and explicitly shown that quark 
mass effects are not suppressed. We point out that it is a feature of the off-forward case 
only and this term is absent in the forward limit. 

We would like to thank A. Harindranath and M. V. Polyakov for helpful discussions. 
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